Binary Black Hole Encounters, Gravitational Bursts and Maximum Final Spin 
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The spin of the final black hole in the coalescence of nonspinning black holes is determined by 
the "residual" orbital angular momentum of the binary. This residual momentum consists of the 
orbital angular momentum that the binary is not able to shed in the process of merging. We study 
the angular momentum radiated, the spin of the final black hole and the gravitational bursts in a 
sequence of equal mass encounters. The initial orbital configurations range from those producing 
an almost direct infall to others leading to numerous orbits before infall, with multiple bursts of 
radiation while merging. Our sequence consists of orbits with fixed impact parameter. What varies 
is the initial linear, or equivalently angular, momentum of the black holes. For this sequence, the 
final black hole of mass Mh gets a maximum spin parameter a/Mh ~ 0.823, with this maximum 
occurring for initial orbital angular momentum L/M\ ~ 1.176 . 
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PACS numbers: 04.60.Kz,04.60.Pp,98.80.Qc 



A few years ago, after a decades-long period of de- 
velopment, breakthroughs were made in computational 
modeling of strong gravitational fields that now allow nu- 
merical relativists to successfully simulate binary black 
holes (BBH) from inspiral through merger. In general 
terms, there are now two computational recipes to fol- 
low. One of them is based on a generalized harmonic 
formulation of the Einstein equations [l|, Q and uses ex- 
cision H || of the black hole (BH) singularities. The 
other recipe, called the moving puncture recipe, involves 
a BSSN [5j, [fj formulation, punctures to model BH singu- 
larities and a gauge condition for these punctures to move 
throughout the computational domain 0, [1] ■ Using these 
recipes, many studies involving interacting BHs and their 
generated gravitational radiation have been carried out, 
including gravitational recoil @, 01, HO, [H, [H GS 0, 
spin hang-up ]1& and matches to post-Newtonian (PN) 
approximations!!?], fisl ]. Most center on astrophysical 
implications and connection to future gravitational wave 
observations. 

BBH simulations also enable studies of strong non- 
linear phenomena regardless of traditional gravitational 
astroph ysic s consequences. A recent example is the work 
in Ref. [19| on the self-similar behavior found in the ap- 
proach to the merger/flyby threshold of BBHs. Similar 
merger thresholds in BBH encounters or scatterings form 
the context for our work. 

We consider orbits in which the BHs initially fly past 
one another, but then fall back to orbit and merge. We 
focus on the gravitational waveform and the angular mo- 
mentum radiated from such encounters. Scrcndipitously, 
we find significant astrophysical implications, both the 
existence of a maximum in the final BH spin and of mul- 
tiple encounter orbits with associated multiple bursts of 
gravitational radiation. Ref. considered only the first 
close encounter or "whirl," and the study did not extend 
the evolutions to find possible fall-back orbits such as 



those here considered. The work in Ref. [19| and our 
work here have to date been the only studies considering 
these highly eccentric orbits; while there have been high- 
order PN studies of inspiral, cases studied so far have 
described relatively smooth inspirals (2fj| . 

All our orbits arc parabolic or hyperbolic encounters. 
Depending on the merger, the fraction of angular momen- 
tum radiated varies significantly (0.05 £ J ra d/L £ 0.55 
with L the initial orbital angular momentum of the bi- 
nary). This emission of angular momentum sets an upper 
limit of a/Mh w 0.823 for the spin parameter of the final 
BH; this maximum occurs when L/M? rj 1.176, with 
Mh the mass of the final merged BH. 

As in our previous BBH studies [111, [2l|, [22j, we use 
a code based on the BSSN formulation and the moving 
puncture recipe. The results here were obtained with 
a 634 3 M computational domain consisting of 10 refine- 
ment levels, with finest resolution of M/52. We set up 
nonspinning equal mass BHs using Bowen-York initial 
data (23|. The mass of each BH is M/2, computed from 
\J Aah/lQn with A a h their apparent horizon area. The 
data have the BHs on the x-axis: BH± is located at ±5 M 
and has linear momentum P± = (=F-P cos#, ±P sin (9, 0). 
We keep the angle constant at 9 = 26.565° = tan (1/2); 
thus the impact parameter is ~ 4.47 M. The total ini- 
tial orbital angular momentum is given by L/M 2 = 
10 (P/M ) sm6z. We obtain a one-parameter family of 
initial data by varying the magnitude of the initial mo- 
mentum in the range 0.1145 < P/M < 0.3093. At the 
lower limit of the momenta, merger occurs within less 
than half an orbit of inspiral. We then consider succes- 
sively higher initial momentum until we find solutions 
that will clearly require a very long, "infinite" time to 
merge. 

The results are summarized in Figs. [T] and [2J The 
top panel in Fig. [1] shows the spin a/Mh of the final 
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BH as a function of the initial orbital angular momen- 
tum L/M%. The spin and mass of the final BH were 
computed using the apparent horizon formula [111, l24l ]. 
The bottom panel of Fig. [T] displays the fraction of an- 
gular momentum radiated (J r ad/L = 1 — aMh/L). Fig- 
ure [2] shows, as a function of L/M 2 , in the top panel 
the final mass Mh/M adm relative to the total ADM mass 
and in the bottom panel the fraction of energy radiated 
E ra d/M adm = 1-M h /M adm . The vertical lines in Figs. Q] 
and [3 denote the value of L/M 2 where a/Mh is maxi- 
mum. We have also calculated both the radiated an- 
gular momentum and energy via the Weyl tensor. The 
results are consistent with those in Figs. Q] and [2] How- 
ever, the values obtained form J ra d/L = 1 —aMh/L and 
E r ad/Madm = 1 — Mh/M a dm are more accurate because 
they are not as susceptible to resolution effects as those 
derived from wave extraction. We have carried out sim- 
ulations at resolutions of M/45, M/48, Af/52 and M/64 
for ten representative cases in Figs. [T] and to check con- 
vergence and make error estimates. We found that the 
results are consistent with the 4th-order accuracy of our 
code and that the errors in the quantities displayed in 
these figures are not larger than 3%. 

We have selected six encounters that are representative 
of the different behaviors in our series. These six cases 
are L/M 2 = 0.512, 1.208, 1.352, 1.376, 1.382 and 1.387 
or equivalently L/M% = 0.521, 1.282, 1.522, 1.480, 1.498 
and 1.554. We will refer to them as encounters Ea, Eb, 
Ec, Ed, Ee and Ef, respectively. Cases Ed, Ee and Ef 
correspond to the last three points in Figs. [1] and [2j 

For L/M% ^,0.8, the radiated angular momentum is 
Jrad/L ;$ 0.15, so the final BH has a/Mh close to L/M 2 . 
The evolution is rather simple in these cases: immedi- 
ate merger, with minimal inspiral. For instance, in case 
Ea (Fig. [3]), L/Ml = 0.521, and J ra d/L = 0.05; thus 
most of the angular momentum goes into the final BH, 
a/Mh = 0.496. Fig. [UEa shows the corresponding ra- 
diated gravitational wave (MrRe^ 2 ). All waveforms 
were extracted at radius 50 M. 

As the initial angular momentum increases, the radi- 
ated angular momentum also increases, suppressing and 
limiting the spin of the final BH. Eventually for large 
enough initial angular momentum, so much angular mo- 
mentum is radiated that, as seen in Fig. [1] the final spin 
reaches a maximum of a/Mh ~ 0.823 at L/M 2 w 1.176. 
Fig. [3lEb shows the tracks of the BHs in the neighbor- 
hood of this maximum. Fig. [¥]Eb shows the correspond- 
ing radiated waveform. For even larger initial angular 
momentum, the spin of the final BH actually decreases 
for increasing L/M 2 . The reason is that the merger is not 
only preceded by several hang- up orbits [H, [l9[ , but also 
the merger yields a highly distorted BH that radiates co- 
piously as it settles down. Case Ec with a/Mh ~ 0.68 and 
L/M 2 w 1.522 represents this situation in which almost 
50% of the initial angular momentum is radiated (see 
path in Fig. OEc and radiated waveform in Fig. [UEc). 

A persistent feature of the mergers with L/M 2 £ 1.3 
is that the separation between the BHs (the coordinate 
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FIG. 1: Top panel, spin of the final BH a/Mh and, bottom 
panel, angular momentum radiated J ra d/L vs the initial or- 
bital angular momentum L/M%. 

distance between the punctures) decreases monotonically 
with time (monotonic inspiral). Comparing cases Ea, 
Eb and Ec in Fig. [4] we see general qualitative agree- 
ment: inspiral-generated gravitational waves with fre- 
quency and amplitude increasing in time, followed by es- 
sentially fixed- frequency ringdown waves. There is, how- 
ever, a hint of disappearance of the monotonic spiral in 
case Ec. The amplitude of the gravitational radiation 
has a "shoulder" at about time ~ 110 M. For a period 
of time equal to two wave oscillations, the decline of the 
amplitude ceases and then recommences. The relative or- 
bital separation as a function of time (Fig. [HlEc) clearly 
shows there is a plateau in the separation centered at 
time ~ 50 M, which is absent for cases Ea or Eb. For 
a brief period of time there is a closely circular phase in 
which the BHs "want" to fly apart, but just manage to 
stay at roughly constant separation. 

The last three points in Figs. [T] and [Hare the cases la- 
beled Ed, Ee and Ef. They describe orbits without imme- 
diate merger but "escape" and recapture; they all show 
initial approaches followed by increasing mid-evolution 
separations of 14 M, 25 M and 42 M before the final 
merger (see Figs. [3] and H]). Because the interaction in- 
volves two close approaches, there are two bursts of grav- 
itational radiation, one from the first flyby [25| and the 
other from the merger (see Fig. [5]). Wc arc currently 
investigating astrophysical implications of detections of 
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FIG. 2: Top panel, mass of the final BH Mh/M a dm and, bot- 
tom panel, energy radiated E ra d/M a dm vs the initial orbital 
angular momentum L/M%- 



FIG. 3: BH tracks of the encounters. The coordinate dimen- 
sions of the top four panels are 12 M x 12 M and 16 M x 16 M 
for the 2 bottom panels. 



these multiple gravitational bursts and hangups in glob- 
ular clusters [261 ] . 

For the Ef case, there is an approximate hangup with 
separation ~ 4 — 5 M around time ~ 950 M similar to 
the shoulder seen in Fig.[H]Ec around time ~ 50 ill. This 
structure shows up in the waveform for this Ef case; we 
actually see a (lower amplitude) precursor to the radi- 
ation burst associated with the merger, a hint that or- 
bits with many repeated bounces are possible. For even 
slightly (0.1%) greater initial angular momentum than 
case Ef, the BHs complete approximately one loop and 
then escape. This is a possible indication of chaotic be- 
havior (exponential dependence on initial conditions, c.f. 
Ref. [lj|). Repeated bounce orbits would have to be 
found with initial angular momentum very slightly above 
that which resulted in Fig. [3]Ef. As with all critical phe- 
nomena, the problem becomes one of careful tuning of 
the parameters. Note that these interactions of nonspin- 
ning BHs produce chaotic orbital dynamics, in contrast 
to the chaos found in spin evolutions [13, H8| . 

One of the main conclusions of our work is that there is 
an upper limit on the Kerr parameter for the final merged 
BH from nonspinning BH merger. For our sequence this 
maximum is a/Mh ~ 0.823 . We can understand this 
observation by examining the timing of the formation of 
the final BH and the radiation from the merger. It ap- 
pears that the merger occurs through an intermediate ex- 
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FIG. 4: Waveforms for the Ea, Eb and Ec encounters. 
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FIG. 5: Waveforms for the Ed, Ee and Ef encounters. 



cited state which is essentially a highly distorted BH. We 
say "essentially" because a substantial amount of angular 
momentum is also radiated in the plunge immediately be- 
fore the apparent horizon forms. This is consistent with 
close limit BBH calculations [29( that show merging BHs 
behaving like a perturbed BH, even before a common ap- 
parent horizon forms, so long as the merging BHs are in- 
side the peak of the effective potential of what will be the 
final BH. This intermediate state emits the largest part 
of the radiated energy and angular momentum. Because 
this mechanism is universal (excitation of such a state 
is inevitable, and it will inevitably radiate), it suggests 
that no merger of equal mass (or presumably, roughly 
equal mass) BHs can lead to a final BH with maximal 
spin parameter a/Mh ps 1. This result does not directly 
affect spin up by accretion since mass accretion will not 
excite the low I modes that strongly radiate angular mo- 
mentum. Thus typical gas accretion can in principle lead 
to final spins much closer to the limit a/Mh = 1. 
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FIG. 6: BBH coordinate separation for the Ec, Ed. Ee and 
Ef encounters. 



